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Abstract 

We report an analysis of the Anderson transition in an SU(2) model with chiral symmetry. Clear single parameter 
scaling behaviour is observed. We estimate the critical exponent for the divergence of the localization length to be 
v = 2.72 ± .02 indicating that the transition belongs to the symplectic universality class. 



Key words: Anderson transition; Spin-orbit coupling; Chiral symmetry 



1. Introduction 

It is thought that the critical phenomena of the An- 
derson transition may be classified according to three 
universality classes (orthogonal, unitary and symplec- 
tic) depending on the symmetries of the Hamiltonian, 
i.e. time reversal and spin rotation symmetries. A dif- 
ferent value of the critical exponent v describing the di- 
vergence of the localization length is expected to char- 
acterize each universality class. 

In the absence of any diagonal disorder the SU(2) 
model Hamiltonian may have chiral symmetry. This 
depends on the boundary conditions and the number 
of sites in each direction. Near the band center chiral 
symmetry affects the localization of electrons, e.g., an 
even-odd system size dependence of the quasi- Id lo- 
calization length is observed [1]. Here we investigate 
whether or not chiral symmetry also affects the critical 
exponent v. 

A system belongs to the symplectic universality class 
if its Hamiltonian commutes with a time reversal op- 
erator T that satisfies T 2 — — 1, e.g., systems with 
significant spin-orbit coupling. Systems in this univer- 
sality class exhibit an Anderson transition even in 2D. 
Recently we analysed the SU(2) model with an on-site 



random potential (diagonal disorder), which breaks 
chiral symmetry, and estimated the critical exponent 
v = 2.73 ± .02 [2]. 

Here we study the SU(2) model without an on-site 
random potential. In spite of the fact that there is no 
random potential, we find that there is an Anderson 
transition at a critical energy E c . Further this energy 
E c is far from the the band center and chiral symmetry 
does not affect the observed critical phenomena. 



2. Model and Method 

The Hamiltonian of the SU(2) model describes non- 
interacting electrons on a simple square lattice with 
nearest neighbour SU(2) random hopping 
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where cJ CT (c; CT ) denotes the creation (annihilation) op- 
erator of an electron at the site i with spin a. We dis- 
tribute hopping matrices randomly and independently 
with uniform probability on the group SU(2). This cor- 
responds to a and 7 uniformly distributed in the range 
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Fig. 1. In A vs Fermi energy E for system sizes 
L — 11, 16, 23, 32, 47, 64. The solid lines arc the fit. 

[0, 2ir), and /3 distributed in the range [0, it/2] accord- 
ing to the probability density, P(/3)d/3 = sin(2/3)d/3. 

We calculate the localization length A on a quasi- Id 
strip whose width is L with the transfer matrix method 
[3] and analyse the dependence of the re-normalized 
localization length A = A/L on the energy and the 
width of the strip [4] . We also calculate the density of 
states with the Green's function iteration method [5] . 

In the course of our calculations we noticed a tech- 
nical problem with the method described in [3]. When 
fitting the numerical data it is important to have not 
only an estimate of the localization length A but also 
an estimate 8 A of how accurately A has been estimated. 
Otherwise a reliable fitting of the data is not possible. 
We have found that the error 8 A is over estimated by 
the method of [3] . To overcome this problem we divide 
the quasi- Id strip into blocks whose lengths arc 10 3 and 
calculate the quasi- Id localization length Ai(i denotes 
the block) in each block and accumulate Ai. From the 
distribution of {A;} we estimate A and 8 A. 



3. Results 

We calculated A as a function of Fermi energy for 
systems with widths L = 11, 16, 23, 32, 47 and 64. The 
accuracy of the localization length data is 0.1%, with 
the exceptions of L = 47, 64 where the accuracy is 
0.2%. To achieve this accuracy systems of length of the 
order of 10 7 to 10 8 are required. We impose the periodic 
boundary conditions in the transverse direction. For 
this choice of boundary condition systems with even L 
have chiral symmetry, while those with odd L do not. 
We fit even and odd data separately first and then in 
combination. The results are displayed in Table 1. The 
estimates of the critical parameters for all three fits 
are in close coincidence and we conclude that chiral 



Table 1 

The number of data Nd and goodness of fit Q. The best 
fit estimates of the critical energy E c , lnA c and the critical 
exponent v with 95% confidence intervals. There arc 7 fitting 
parameters. 





N d 


Q 


E a 


In A e 


V 


Even 


159 


0.1 


3.2531±.0002 


0.613±.002 


2.72±.02 


Odd 


164 


0.2 


3.2531±.0002 


0.613±.002 


2.72±.03 


Both 


323 


0.1 


3.2531±.0001 


0.613±.001 


2.72±.02 
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Fig. 2. The DoS near the band edge for a system of 
size 100 X 100. The vertical line shows the mobility edge 
(E c ~ 3.2531). (The imaginary part of the energy is 0.005.) 

symmetry does not change the universality class. In 
Fig. 1 we show the result for the combined even and odd 
data. The lines cross at single point, indicating both 
that corrections to scaling due to irrelevant variables 
are negligible and also that an even-odd effect does not 
exist. The estimates v and lnA c are consistent with 
those in [2] for the symplectic universality class. 

The density of states for a system of the size 100 x 100 
averaged over 10 3 ensembles is shown in Fig. 2. In the 
critical region, the density of states changes slowly. 
This behaviour is different from that of the random 
magnetic flux U(l) model in 3D where the density of 
states changes rapidly in the critical region [6]. This 
might be the reason why we obtain clearer single pa- 
rameter scaling here than in a study of the critical phe- 
nomena of the U(l) model [7]. 
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